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ABSTRACT 



Aims. The purpose of this paper is to characterize the statistical properties of solar granulation in the photosphere and 
low chromosphere up to 650 km. 

Methods. We use velocity and intensity variations obtained at different atmospheric heights from observations in Ball 
4554 A. The observations were done during good seeing conditions at the VTT at the Observatorio del Teide on Tenerife. 
The line core forms rather high in the atmosphere and allows granulation properties to be studied at heights that have 
been not accessed before in similar studies. In addition, we analyze the synthetic profiles of the Ban 4554 A line by 
the same method computed taking NLTE effects into account in the 3D hydrodynamical model atmosphere. 
Results. We suggest a 16-column model of solar granulation depending on the direction of motion and on the intensity 
contrast measured in the continuum and in the uppermost layer. We calculate the heights of intensity contrast sign 
reversal and velocity sign reversal. We show that both parameters depend strongly on the granulation velocity and 
intensity at the bottom photosphere. The larger the two parameters, the higher the reversal takes place in the atmo- 
sphere. On average, this happens at about 200-300 km. We suggest that this number also depends on the line depth of 
the spectral line used in observations. Despite the intensity and velocity reversal, about 40% of the column structure 
of granulation is preserved up to heights around 650 km. 
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1. Introduction 
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Convective motions in the solar atmosphere have been stud- 
ied in detail during the past decades, using different types of 
observations, their analysis, and interpretation. The sum- 
mary of the in v estiga tions before 2 000 c an be found in 
lEspagnet et"all (|1995D and iBrandtl (|2000f) . According to 
the final results, the photosphere can be divided into two 
regions with different physical properties with the border 
between them lying at heights about 170 km above the 
continuum formation height (see lEspagnet et al.l . I1995L for 
details). This conclusion was reached based on the ampli- 
tudes of intensity and velocity fluctuations first decreasing 
with height, reaching minimum at about 170 km and then 
starting to increase again. The intensities and velocities at 
heights above and below 170 km do not se em to be corre- 
lated with each other. On the other hand, lEspagnet et al.l 
(|199 5^ came to the conclusion, based on their own inves- 
tigations, that the photosphere can be divided into these 
two regions only in intensity fluctuations, since the vertical 
convective velocity structure is preserved until heights as 
high as 500 km. A s imilar result has also been obtained by 
ISalucci et all (|1994f ). 

During the past few years significant efforts have been 
made to determine the precise height where the inversion of 
intensity contrast (or temperature) takes place. The results 
by different authors diverge significantly in this respect. 



Values as low as 60 km are found by iKneer et all (|1980f) 
from spectrogram s of the Mg i b2 line ta.ken a t the quiet sun 
disk center. While iBendlin fc Vol kmeil (119931 ) give values as 
high as 270 km from 2D spectral images obtained with tun- 
ing a Fabry-Perot interferometer in the Fe i 6303 A line in a 
moderately active region near disk center. Different values 
are also given in the more recent liter ature, e.g. 250 km 
fe ostyk fc Sh chukina. '20 04D, 20 km (Janssen & Gauzzi 
200i), 170 km Puschman net all ([2005). In the last paper, 
the authors perform a careful analysis and NLTE inversion 
of the observed spectral profiles of the several photospheric 
Fc I lines and present the resulting temperature and veloc- 
ity stratifications as functions of both optical depth and 
geometrical height. In addition to the height of inversion 
of the intensity contrast (temperature), they also report 
that only the large-scale structures (about 4") penetrate 
to higher layers without the lost of correlation though the 
whole atmosphere. Sim ilar results were also reported by 
iPuschmann et all (|2003i) . where the authors find that only 
the intensity structures with sizes over 2" at heights about 
435 km are still connected with those at the continuum 
level. 

In all the p apers reported above (except 
iKostvk fc Shchukinal . |2004[ ). the vertical velocities were 
observed in phase though the whole photosphere pen- 
etrating into the highes t layers under study. Only 
iKostvk fc Shchukinal (|2004[ ) have obtained evidence of the 
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velocity reversal at heights above 490 km for the elements 
with the highest contrast in their data (above 6-9%). 

The spectral lines used in previous studies gave the pos- 
sibility of analyzing statistical properties of solar granu- 
lation mainly in photospheric layers below the tempera- 
ture minimum. The aim of the present paper is to extend 
this analysis to the higher layers taking advantage of spec- 
tral observation s in th e B an 4554 A line. Acc ording to 
lOlshevskv et al.l (|2008D and lSiitterhn et al.l (|200lD . the core 
of this line forms rather high at about 700 km. It makes the 
line particularly useful for studying the vertical structure of 
the solar atmosphere, since we can follow this structure to 
much greater heights tha n in the previous analys is. Recent 
theoretical investigation ijShchukina et al.l l2009l hereafter 
referred to as Paper I) based on NLTE modeling the Ban 
4554 A line in a three-dimensional (3D) hydrodynamical 
model have shown that this line is a valuable tool for the 
Doppler diagnostics along the whole photosphere and lower 
chromosphere. 

In the present study we use theoretical understanding 
attained in Paper I for interpreting the new high-resolution 
spectral observations in the Ban 4554 A line. We investi- 
gate the column structure of granulation up to the heights 
above the solar temperature minimum, derive the height of 
the intensity contrast sign reversal and the change of sign 
of the vertical motions, and propose a 16-column model of 
granulation, depending on the direction of motion and in- 
tensity in the low and high photospheres. The results of 
observations are compared with the NLTE calculations of 
the Ban 4554 A line profiles in the 3D HD model of solar 
convection bv lAsplund et al.l (|2000D . 



2. Observations 

Observations were obtained in July 21, 2004 at the 
70-cm German Va cuum Tower Telescope (VTT, see 
ISchroter et al.l . Il985l . for description) at the Observatorio 
del Teide in Tenerife during the good seeing conditions. 
The observations were done with Echelle spectrograph. The 
spectral region was centered on the Ban 4554 A line. Using 
the narrow-band Can K and Ha slit-jaw images, we se- 
lected a quiet region free of magnetic activity close to the so- 
lar disk center. Spectral images were recorded with a CCD 
camera with 1024 by 1024 pixels. The size of the images in 
spectral direction was 4.1 A. The width of the spectrograph 
entrance slit was 100 /im (equivalent to 0."46). One pixel 
corresponds to 0."087 on the solar surface. This way, the 
field of view of the telescope was 0."46 by 89". However, 
after the reduction, only about 730 pixels were found to 
be useful for the analysis, equivalent to the 63. "5 of the 
solar surface. The time series of Ban 4554 spectral images 
was obtained with a fixed slit position of the total duration 
of 69.9 minutes with an interval of 7 seconds. The expo- 
sure time was 1 second, and a total of 600 exposures were 
taken. The image of the observed area was stabilized using 
the adaptive optics (AO) system guiding on solar granula- 
tion. The average Fried's parameter, as given by AO, was 
fluctuating around Rq = 7 — 8. The spatial resolution of ob- 
servations limited by the image motion on the spectrograph 
slit due to earth atmospheric turbulence was estimated to 
be about 0."3-0."7. 



3. Reduction of observations and calculation of line 
parameters 

Following the standard proced ure, all 600 s pect ral images 
were reduced for flatfield {e.g. iKiselmanl . Il994i ). The flat 
field images were recorded immediately after the time se- 
ries of Ban spectra in the same wavelength at the solar 
disk center. The spatial structure of the solar surface was 
removed from the flatfield images by rapid motions of the 
telescope scan mirror. The reduction for flatfield included 
the following steps: 

— removing the dark current out of the spectral images 
and the flatfield images, 

— correction of all the images for the inclination and dis- 
tortion of the spectrograph sit with respect to the pixel 
direction of the CCD, 

— removing the solar spectral lines from flatfield images, 

— splitting the fiatfield images into the two components. 
Here we need additional explanations. Both CCD cam- 
era and the spectrograph slit have some dust on them. 
The dust at the CCD camera appears as dark dots on 
the spectral images, while the dust at the slit appears 
as dark horizontal lines, similar to intergranular lanes. 
As the dust position on the slit can change with time, 
the position of these dark lines on the spectral images 
and on the flatfield images is not always the same. Thus, 
we had to split the flatfield images into the images con- 
taining only a CCD dust component (matrix fiatfield) 
and those containing only the slit dust component (slit 
fiatfield). 

— removing the matrix flatfield from spectral images. This 
can be done automatically. Removing the slit fiatfield re- 
quires individual control of all the spectral images, since 
the positions of the dark lines do not always coincide. In 
some cases we had to recalculate the slit fiatfield images 
by shifting them in a spatial direction by a fraction of 
a pixel. 

The next step in the data reduction included calculat- 
ing of the intensity and velocity variations at different po- 
sitions of the Ban line profiles. For that we used a A-meter 
te chnique. The description of thi s method can be found 
in IStebbins fc Goodd (|1987D and IShchukina et al.l (|2009l . 
Paper I). 

We used 55 reference widths of the average Ban 4554 
A line intensity profile AAiv = A A,, — AAf,, where AA^ 
and AAfc are the wavelengths of the red and blue wings, 
respectively. For each individual Ban profile at each spa- 
tial position X (varying from 1 to 730) and time moment 
t (varying from 1 to 600), we then found intensity vari- 
ations dI{t,x,W) and velocities V(t,x,W) at the corre- 
sponding 55 reference widths. We slightly modified the 
standard A-meter procedure as described below, keeping 
in mind that the formation heights of the blue and red 
wing intensity points of the line profile are, in general, not 
the same becau s e of the Doppler velocity shift (see e.g. 
IShchukina et al.l . l2009l . Paper I). The velocities were calcu- 
lated separately in the red and blue wings of the Ban line 
in the following way: 

Vr{t,X,W) = SVr{t,X,W)+V{W) (1) 

Vb{t,x,W) = dVb{t,x,W)+V{W) 
where 

dVrit,x,W) = (c/A) • SXr (2) 
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Fig. 2. Convective (left panel) and oscillatory (right panel) components of the velocity field corresponding to heights 
about 650 km as a function of the slit position and time. 



300 




Velocity, km/s 

Fig. 1. Bisector due to hyperfi ne and isotop i c shif t calcu- 
lated in the MACKKL model fMaltb v et all . |1986[ ) (small 
open circles), and bisector of the observed profile, average 
over space and time, caused by the nonthermal velocity field 
(solid hue). 

6Vbit,x,W) = {c/\)-S\b. 

In these equations, F(VF) is the bisector of average profile 
due to to nonthermal velocity field 

V{W) = (c/A) • (AA; + Ah)/2 - VuFSi , (3) 

where Vhfsi is the line bisector caused by the hyperfine 
structure (HFS) and isotopic shift, c velocity of light and, 
A the line wavelength. Furthermore, in eq. ((2]) SXb^r = 
AXb^r — ^^b,r denotes the wavelength shift of the blue (red) 
intensity point belonging to a certain spectral width rela- 
tive to the corresponding blue (red) point of the average 
profile. 

We derived the bisector Vhfsi using the Ban 4554 A 
line profile calculated in the MACKKL model (see e.g. 
IShchukina et al.L l2009l Paper I). Figure [1] shows that the 



bisector T^psi is negative (between —0.3 and —0.2 km/s), 
while the nonthermal velocity bisector y(VF) of the aver- 
age observed profile changes sign ranging from —0.1 km/s to 
« 0.25 km/s. The velocity values were corrected for the dis- 
placement due to the Earth rotation and revolution. Note 
that through this paper negative velocities mean downflows 
and positive velocities mean upflows. The only exception is 
Fig. [T] where, in agreement with Fig. 2 of Paper I, we as- 
sume negative velocity direction to be toward the observer, 
corresponding to an upflow. 

In a similar way, we calculated intensity variations (con- 
trast) dl{t, X, W) relative to the intensity I{W) of the Ban 
line profile averaged over space and time: 

6I{t,x,W) = iI{t,x,W) - I{W))/I{W) . (4) 

A through analy s is of t he Doppler diagnostic techniques 
bv IShchukina et all (|2009l Paper I) has shown that the Vb 
and Vr velocities calculated from the line profile of the 
Ban 4554 A line are, in fact, very close to each other. 
According to their analysis, both Vb and Vr correspond bet- 
ter to heights of formation Hb of the blue wing intensities, 
rather than to Hr of the red wing intensities. Guided by the 
results of Paper I we use the average values V = {Vb + Vr)/2 
throughout the present paper and ascribe their origin to the 
average heights of formation of the blue wing intensities. 

Variations of intensity SI{t,x,W) and velocities 
Vr{t, X, W) and Vb{t, x, W) are mainly caused by oscillatory 
and convective motions. By convective motions we mean 
the motions of individual granules and intergranular lanes. 
To separate oscillatory and convective components we per- 
formed filtering in the Fourier domain based on the diagnos- 
tic fc — a; diagram. According to this diagram, we limited 
the wave motions in temporal frequency between to =1.8 
and 5.7 mHz and convective motions by frequencies below 
iu = 2.2 mHz. Additional filtering was performed for con- 
vective motions in the spatial frequency, leaving only varia- 
tions wi th k less than 0.18 Mm~ ^. Si milar filtering was ap- 
plied in iKhomenko et al.l (|2001l ) and iKostik fc Khomenkcl 
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Fig. 3. Vertical velocities along one of the slices in the 
3D model. Negative velocities correspond to intergranules 
(downflow, blue color) , while positive velocities correspond 
to granules (upflow, red color). 



([2003)- As the frequency domains overlap, the wave mo- 
tions can, in principle, contribute to convective motions. 
The maximum contribution of the wave motions at w = 3.3 
mHz reaches 4.2%. We checked that decreasing this con- 
tribution to 3.1% (more narrow filter) or increasing it to 
20% (more wide filter) produces only an insignificant in- 
fluence on our results. As an example we show in Fig. [2] 
the convective (left) and oscillatory (right) components of 
the velocity corresponding to heights around 650 km as a 
function of space and time. The vertical axis represents the 
velocity amplitudes. The different temporal and spatial be- 
havior of the both components is evident from the figure. In 
the rest of the paper we focus our analysis on the convective 
component of the velocity and intensity variations. In addi- 
tion, we only select a part of the time series (of about 200 
exposures) when the seeing conditions were the best; i.e., 
the estimated spatial resolution was not worse than 0."4. 

4. Formation heights along the observed Ba ii line 
profiles 

The formation heights of the 55 reference widths of the 
Ball 4554 A intensity profiles were obtained from NLTE 
calculations in a 3D hydrodynamic model atmosphere by 
(jAsplund et al.l . r200Cl[ ). Detailed description of the Barium 
atomic model, NLTE calculation s, and the 3 D hyd rody- 
namical snapshot is g iven in Ol shevskv et al.l ()2008[ ) and 
IShchukina et al.l ()2009l Paper I) . In this section we present 
additional information that can be useful for understanding 
our results. 

Figure [3] visualizes vertical variations in the velocity 
along one of the slices of this snapshot. The geometrical 
height is the same in all the points of the snapshot, but the 
optical depth scale is fluctuating from column to column 
due to horizontal inhomogeneities of the thermodynamic 
quantities. The zero point of geometrical height H = is 
assigned by averaging over all grid points of the snapshot 
the locations where the c ontinuum optical depth t^qoo ~ 1 
(see lAsplund et al] . l2000f) . 

The vertical extent of the snapshot was optimized to 
reach high enough layers (till 900 km) to allow reliable 
calculation of the cores of strong lines (see lAsplund et al.l . 



[2000l:[Stein"fc NordlundfTOOStlNordlund fc Dravinl [iQQOl ) . 

Taking into account that the core of the Ban 4554 A line 
i s formed more than 1 00 km below the upper boundary 
(jShchukina et al.l . l2009l Paper I), one can expect the ef- 
fects of the upper boundary conditions to be negligible for 
spectral synthesis of this line. An excellent agreement of 
the spatially average synt hetic profile and observed profile 
from the Liege atlas (Delbouille et al.Lll973D . as well as pro- 
file obtained from obser vations at the VTT, su pports this 
conclusion (see Fig. 4 in IShchukina et all . l2009l Paper I). 

The solar disk-center profiles of the Ban 4554 A line 
were calculated by means of ID multilevel NLTE radiative 
transfer employing a set ID models from individual grid 
points of the 3D snapshot as input. Such an approach is 
known in the literature as 1.5D approximation. 

For each grid point (ID model), we calculated forma- 
tion heights of the blue wing intensities Hi, corresponding 
to 55 reference widths AXyy of the Ban profile. Following 
Shchukina et al. (20091, Paper I), we used the concept of 
Eddington-Barbier heights of line formation; i.e., we evalu- 
ated heights where the line optical depth at a given wave- 
length point AA is equal to unity: t(AA) = 1. We then 
assumed that the information on vertical velocity and inten- 
sity at wavelength point AA comes from this height. We re- 
peated the calculations of Hi, heights for the all grid points 
of the 3D snapshot and then computed an average value 
(H) for each of the 55 reference widths AAi^. These forma- 
tion heights were ascribed to the corresponding widths of 
the observed intensity profiles. Thus, the formation heights 
calculated this way are model-dependen t. 

According to'Shch ukina et al] (|2009l Paper I), the mean 
formation heights {H) of each section of the Ban 4554 A 
line profile have a well-pronounced dependence on the spec- 
tral width of this section. Thus, the intensity and velocity 
information coming from each width level AAvf of the ob- 
served profile roughly corresponds to a different height in 
the atmosphere. The selected widths AXw of observed pro- 
files vary between 280 mA and 30 mA. The corresponding 
range of the mean formation heights lies between — 25 km 
and ~ 650 km. 

The lowest height level corresponds to the width 
AXw = 280 mA, and its average formation height is 
(H) = — 25 km. We used the intensity variations at this 
level as a criterion for separating granular and intergranular 
regions. Taking into account that this height is expected to 
be close to the formation height of the continuum intensity, 
we call "granules" those features whose intensity I{t, x, W) 
at AXw _— 280 mA {{H) — —25 km) was larger than the 
average I{W) . The classical definition of "granule" and "in- 
tergranule" can only be applied to the continuum formation 
heights. Higher in the atmosphere, this definition does not 
make sense. 



4.1. Spatial smearing and determination of Fried 's parameter 

To directly compare between the synthetic and observed 
spectra, we performed a convolution of the synthetic pro- 
files with a corresponding modulation transfer function 
(MTF) representing the action of the 70-cm telescope 
and Earth's atm osphere (see e.g. iKhomenko et al.l . 120051 : 
IShchukina et al.l . [2009) . In such an interpretation the qual- 
ity of the resulting image depend s on the tele scope diame- 
ter D and Fried's parameter Rq (jFriedl . [19661 ) . To perform 
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Fig. 4. Root-mean-square continuum contrast (5/r.ni.s. as 
a function of Fried's parameter Rq for wavelength A = 
4554 A and for telescope diameter D — 70 cm. The 
Slr.m.s. curve was calculated using 3D radiative trans- 
fer formal solution for cont inuum intensity obtained by 
iTrujillo Bueno fc Shchukinal (|2009f ) in the 3D HD model. 
Filled circles correspond Rq = 8 cm {6Ii ,„ s = 6.1 %) and 
Ro = 38 cm {SI,,n,.s. = 16.5 %). 



the smearing, the original 2D intensity maps at each indi- 
vidual wavelength were Fourier-transformed and multiplied 
by MTF calculated for known telescope diameter (70 cm, 
VTT) and Fr ied's parameter Rp. The corresponding MTFs 
are given m IShchukina et all (|2009[) . An inverse Fourier 
transform was performed to obtain the smeared images, 
i.e., the images affected by the diffraction on the telescope 
ap erture and th e seein g effects. 

iRicort et al.l (|1981f ) proposed determining Rq from the 
observed r.m.s. continuum c ontrast of solar granul ation 
^-^r.m.s. Following their idea, IShchukina et al.l ( 20091 ) cal- 
culated theoretical dependences between ^/r.m.s. and the 
Fried's parameter Rq at different wavelengths. These au- 
thors used 3D radiative trans fer solution for the solar con- 
tinuum intensity obtained by iTruiillo Bueno Sz Shchukinal 
(j2009r ) in the same snapshot of the 3D HD model. We 
applied their theoretical calibration curve at wavelength 
A = 4554 A to define Fried's parameter Rq from our obser- 
vations (see Fig. Sj). We found that, with the spatial res- 
olution reached during the observations (~0."7— 0."3), the 
root-mean-square continuum contrast (5/r.m.s. at the wave- 
length of the Ban line is between 5.4% and 6.3%. Thus, we 
estimate Fried's parameter Rq in the range between 7 cm 
and 8 cm. 



5. Results of observations and comparison with 
model calculations 



According to IShchukina et ah! (j2009l Paper I), the A-meter 
technique applied to the Ball 4554 A line profiles recov- 
ers velocities at different heights rather well. Nevertheless, 
there are three important points that we have to keep in 
mind. 

First, under seeing conditions corresponding to our ob- 
servations, (i?o = 7 — 8 cm) the velocities derived from the 
inner wings (AAvk < 150 mA) of the Ban 4554 A line pro- 
files are more reliable than from the outer wings (see Paper 
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Fig. 5. Top: fraction of the spatial locations where the in- 
tensity contrast sign (squares) and the velocity direction 
(triangles) at a given height (or AA^i/) changes along the 
slit. Bottom: same, but averaged over the velocity and in- 
tensity events. 



I, Fig. 9). In the former case, the correlation coefficient 
between the original and A-meter velocities varies between 
0.7 — 0.8, while in the latter it becomes lower (about 0.6). 

Second, such a parameter as the average formation 
height should be used with caution. In reality, the inten- 
sities corresponding to a single value of AAh^ arc formed 
in a wide range of heights, depending on the model atmo- 
sphere. The peak-to-peak formation height variations can 
be as much as 200 km in the outer blue wings and even 
more (~ 500 km) in the inner blue wings (see Paper I, Fig. 
13). 

Third, the absolute values of the A-meter velocities are 
lower than the original ones. This is not surprising. On 
the one hand, we recover velocities from the synthetic pro- 
files degraded by atmospheric seeing. This should lower the 
absolute values of the velocities. On the other hand, the 
A-meter technique itself leads to smoothing and vertical 
averaging of the velocities within the heights where the line 
intensity is formed. 

Taking these limitations into account, we discuss below 
in this section the result s of application of A-meter tech- 
nique to our observations. iKostvk fc Shchukinal (|2004f ) have 
shown that the granular motions in the lower photosphere 
preserve their column structure up to height of 490 km. 
We can verify this result because the Ban line allows us 
to investigate the properties of the granulation at greater 
heights up to 650 km. Using 200 spectral images taken dur- 
ing the best seeing conditions of the series, we calculated 
the number of "convective columns" at each of 55 AAw po- 
sitions. This number was defined as a number of contrast 
sign changes (or, independently, the velocity sign changes) 
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Fig. 6. The fraction of the convective structures with a 
given size at different heights in the atmosphere. The curves 
(from the upper to the lower) correspond to average heights 
\h) = -25, 60, 300, 450, and 650 km. Top panel: size cal- 
culated from the velocity variations; bottom panel: size cal- 
culated from the intensity variations. 

along the slit at each height. We normalize the number of 
these columns (independently for velocity and intensity) to 
the total number at the lowest layer. The results are pre- 
sented in Figini The horizontal axis of the both panels of 
this figure gives the average NLTE formation heights cor- 
responding to 55 values of the profile widths AXw- The 
vertical axis of the top panel gives the percentage of the 
locations of intensity contrast (squares) and the velocity 
direction (triangles) changes along the surface, calculated 
with respect to the lowest layer. The bottom panel of the 
figure gives similar information but averaged over the ve- 
locity and intensity. As expected, the number of the loca- 
tions where the contrast sign or the direction of motion 
changes decreases with height, both for velocity and inten- 
sity events (top panel of Fig. [5]) . There are about 20% fewer 
events for intensity at all heights. However, an unexpected 
result (at least for us) was that about a half of the con- 
vective structures (40%) registered in the continuum are 
maintained with height as high as at 650 km. 

Next we investigated how the size of the convective 
structures changes with height. We define the size as a sim- 
ple difference in spatial coordinate between the successive 
changes in the sign of the contrast or in sign of the velocity. 
The results are shown in Fig. [51 The horizontal axis of the 
figure gives the size of the convective structures, and the 
vertical axis gives the fraction of the structures of a given 
size in percent. The curves are normalized to the corre- 
sponding number of elements found at the lowest layer from 
the velocity variations. The upper curves in both panels in 
Fig. [HI correspond to the average height of (iJ) = —25 km 



and the lower curves to the height of 650 km. The maxi- 
mum of both velocity and intensity distributions appears at 
650-700 km granular size. The number of structures with 
this size decreases with height from about 11% at —25 km 
to ^ 4% at 650 km (as defined by velocity variations) and 
from ~ 9% to ^ 3% (as defined by the intensity varia- 
tions). Convective structures that have larger sizes (about 
1500-1700 km) at the height of continuum formation almost 
always reach the upper layers. Their fraction with height is 
almost constant. It means that the average size of convec- 
tive elements gets larger with height. Here we need to give 
additional explanations. We did not follow the behavior of 
each individual structure with height. It means that, most 
probably, not all the large-size convective structures with 
sizes above 1500 km at (H) = —25 km reach a height of 
650 km. It can instead be a result of merging of the con- 
vective structures with smaller sizes. In either case, we can 
conclude that, in the height range 0-650 km in the solar 
atmosphere, the fraction of the convective structures with 
sizes 1500-2500 km is almost constant. 

How do individual properties of the convective elements 
change with height? We assume that the convective ele- 
ments observed at the uppermost layer at 650 km can ei- 
ther preserve their contrast and velocity direction from the 
continuum formation height or change them by the oppo- 
site. This way we define 16 types of motion. These types 
of motion are outlined in the lower panel of Fig. [71 The 
sign stands for convective elements moving upwards 
or those whose contrast is above the average. The "— " sign 
stands for those elements moving downwards or having the 
contrast below the average. All 16 types of motions are 
indeed present in the solar atmosphere. The upper panel 
of Fig. [3 gives the number of cases corresponding to each 
of the 16 types of motion between and 650 km heights, 
in percent (solid line). The number of cases changes from 
the lowest value of 1.4% (case number 6; downward mov- 
ing hot material in the continuum and upward-moving cold 
material at 650 km) to the highest value of 15.1% (case 
number 14; upward-moving hot material in the continuum 
and upward-moving cold material at 650 km). Only 21.4% 
of the elements show the pure convective character of mo- 
tion {i.e. hot material rises, cold material sinks). The latter 
number contains the sum of the cases 1 and 2 (8.3 -I- 13.1). 
This number can be even lower if the convective motions 
change their direction of motion or sign of contrast several 
times between and 650 km. In contrast, in some 9.9% 
of the cases, the cold material rises and the hot material 
sinks (cases 3 and 4). The remaining 68.7% of the convec- 
tive elements either change their sign of contrast (31.7%) 
or direction of motion (18.8%) or both (18.2%) between 
and 650 km. 

The middle panel of Fig. [7] gives the properties of the 
convective motions in the low photosphere close to the 
continuum formation level (solid line) between and 60 
km. Most of the convective elements preserve their prop- 
erties within this low height range. However, even in this 
case only 75% of the elements follow the classical convec- 
tion behavior, while about 15% of the relatively cold mate- 
rial moves upwards and 10% of the relatively hot material 
moves downwards. 

Similar analysis was also performed with synthetic line 
profiles (see dotted curves in the upper and middle panels 
of Fig. [7t. It is ev i dent t hat the 3D hydrodynamical model 
of lAsplund et al.l (|2000D describes many properties of the 
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Fig. 7. Convective motions in the solar atmosphere accord- 
ing to 16-column model. Upper panel shows the number of 
cases corresponding to each of the 16 types of motions (as 
defined at the bottom panel) between and 650 km. Middle 
panel gives the same but for heights between and 60 km. 
Solid line: observations. Dotted line: 3D HD model. 

observed variations in intensity and velocity, both qualita- 
tively and quantitatively. Only the case of motion number 
10 has a highest percentage in the simulations (15%) at 
heights 0-60 km compared to its absence in the observa- 
tions (middle panel). 

The results shown in Fig. [7] are unlikely to be produced 
by noise in the observations or the influence of the up- 
per boundary conditions in the simulations. By applying 
the filter to separate the wave and the convective motions, 
the high-frequency noise component is completely excluded 
from the observed variations. As for the simulations, the 
Ball line is formed about 100 below the upper boundary of 
the simulation domain, so its influence should be negligible. 

According to the conclusions derived from Fig. [71 the 
convective elements not only change their sign of contrast 
with height but also change the direction of the velocity. 
The latter property has not received much attention in the 
literature, except the work of Kostyk & Shchukina (2004). 
What are the typical heights where the sign reversal of the 
convective intensity and velocity takes place? From the first 
look at our observations, we find that there is a big scatter 
in such heights, from 50 to 650 km. However, we find that 
the cause of such a big scatter is not in systematic or ran- 
dom errors. This scatter comes from the intrinsic proper- 
ties of convective elements, i.e. their contrast and velocity. 



This conclusion is illustrated in Figs. [5] and [H for observa- 
tions with the best r.m.s. continuum contrast SIr.m.s. ~ 6 % 
(solid curves with filled circles). The solid curve in the left 
panel of Fig. [5] shows that the height of the velocity sign 
reversal depends on the strength of the granular veloci- 
ties at the lowest leveQ- The higher is the velocity of gran- 
ules, the higher in the atmosphere the velocity sign reversal 
takes place. The curve is asymmetric with respect to zero 
velocity value. Given the same absolute value of the veloc- 
ity, the downward-moving convective elements change their 
sign of velocity at substantially lower heights than the up- 
ward moving elements. The average value of the height of 
the velocity sign reversal is 212±211 km. 

The right panel of Fig. [8] gives the velocity sign reversal 
height as a function of contrast corresponding the most bot- 
tom layer {{H) = —25 km). Similar to the above, the higher 
the contrast of granules in the continuum, the greater the 
height of the velocity sign reversal. There is again an asym- 
metry in the granulation properties. The bright granules 
change their sign of velocity at much greater heights than 
the dark intergranular lanes. The difference can be as much 
as 300 km for the features with the continuum contrast of 
±7.5%. 

Figure [H illustrates the dependence of the height of con- 
trast sign reversal on the velocity (left panel) and intensity 
(right panel) at (H) — —25 km. The left panel of this fig- 
ure shows that, the lower the downward convective veloc- 
ity at bottom level, the greater the height of the contrast 
sign reversal of these elements. This is what expected con- 
sidering that the observed correlation between convective 
velocity and int ensity in the deep photosp heric layers is 
rather high (e.g. iKostvk fc Shchukinal . I2004D . On average, 
the contrast sign reversal takes place at (H) = 329 ±50 km, 
which is higher than the velocity sign reversal by ~ 100 km. 
According to the right panel of Fig. [51 the contrast sign re- 
versal occurs approximately at the same height both for 
granules and intergranules, as defined by their continuum 
intensity. The higher the absolute value of the continuum 
intensity, the higher the reversal occurs in the atmosphere. 

In general, the agreement between the observed depen- 
dences and the synthetic ones calculated for the same seeing 
{Rq = 8 cm) is fairy satisfactory. However, the agreement 
is worse for results presented in the left panel of Fig. [31 
where the contrast sign reversal height of the elements with 
downward granular velocities calculated from 3D model are 
considerably lower. The reason for that lies in the poor cor- 
relation between the velocity and contrast at the bottom 
layer (H) = —25 km of the 3D model. In the 3D model, the 
correlation coefficient is rather low (~ 0.35), while it is ap- 
preciably higher, (~ 0.65) in observations. These values of 
the correlation coefficients agree with the re sults obtained 
earlier from three iron lines bv Kostvk &: Shchukina (2004). 

Figure [10] gives additional illustration of the granular 
velocity and intensity behavior with height. It shows sev- 
eral examples of the observed horizontal variations of the 
contrast (left) and velocity (right) along the solar surface at 
different heights. The curves demonstrate clearly that the 
amplitudes of the variations becomes progressively weaker 



^ As discussed in Sect. 4, the heights are defined with respect 
to the average location of rsooo — 1 over all grid points of the 
3D snapshot. The velocities and intensities at (H) — —25, as 
given on the horizontal axis of Figs. 8 and 9, are measured from 
the observed profiles at AXw ~ 280. 
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with height and the sign reversal takes place. Such a rever- 
sal is especially evident at spatial locations around 1.5, 3, 
and 5 Mm. 

The dependences discussed above from Figs. [SHTOlwere 
obtained from the ground-based observations under good 
seeing conditions. To what extent are such dependences 
sensitive to the degrading effects caused by Earth's atmo- 
spheric turbulence? To answer this question, we applied 
the same analysis to the set of synthetic profile smeared 
with two values of the Fried's parameters Rq. The results 
for i?o = 8 cm (resolution similar to our observations) are 
shown by open circles with dashed curve in Figs. [5] and [51 
The results for Ro — 38 cm (excellent spatial resolution 
rarely reached in ground-based observations) are shown by 
crosses and dotted curve in the same figures. We considered 
also the case of the profiles with the numerical spatial reso- 
lution of the 3D model when the smearing due to seeing and 
instrumental effects is absent. The latter case is displayed 
in Figs. [TT] and [T^ (dash-dotted curves with open circles). 

The results obtained with the synthetic profiles 
clearly demonstrate that the observed dependences are 
rather sensitive to th e seeing conditions. According to 
IShchukina et al.l (|2009l . Paper I), the spatial smearing re- 
duces the absolute values of the A-meter velocities, partic- 



ularly of the strong downward ones. This explains the dif- 
ference in the continuum velocity scales on the horizontal 
axis of Figs. [21 [HI HH andfT^ (left panels) between the cases 
with i?o = 8, i?o = 38, and the original numerical resolu- 
tion. In the case with Rq = 38 and in the case of the orig- 
inal resolution, the curves have rather symmetric behavior 
with respect to the zero velocity point with nearly twice 
smaller number of the reversal events above the downflow- 
ing areas. Under the resolution similar to our observations 
(i?o = 8 cm), the information on the velocity sign reversal 
above the downflowing areas is lost. The contrast of gran- 
ulation also decreases with decreasing the parameter Rq. 
As a result, all dependences on the granulation contrast 
become smoother and less pronounced under better seeing 
conditions. As follows from Figs. [51 [HI and Figs.[TTl[T3 the 
average height of the velocity sign reversal increases with 
improving spatial resolution, whereas the average height of 
the contrast sign reversal becomes smaller. The difference 
can reach about 100 km between the case with i?o = 8 cm 
and the case with original numerical spatial resolution. 

The evidence presented above gives us confidence that 
our results on the velocity and contrast sign reversal are not 
artifacts caused by specific analysis procedure. Additional 
evidence is presented in Figs. [Til [121 (solid curves with 
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Fig. 10. Left: Observed variations of the granulation contrast along the solar surface at several heights. The curves 
numbered from 1 to 6 correspond to the average heights: (H) = 0, 60, 300, 450, 600, and 650 km. Right: same for the 
velocities. 




Fig. 11. Height of the velocity sign reversal as a function of the velocity (left) and as a function of contrast or temperature 
variations 5T (right) aX H — —25 km. Thick solid lines with filled circles result from 3D model, and dashed lines with 
open circles result from synthetic profiles for the case of no smearing. In the case of the synthetic profiles, the velocities 
and contrasts given on the horizontal axes are the values measured from the profiles at AAi^ = 280 mA corresponding 
to the average height {H) = — 25 km. In the case of the 3D model, we take the velocities and the temperature variations 
at the single layer H = —25 km of each model grid point. 



filled circles) . These data illustrate the velocity sign reversal 
heights and the temperature variation sign reversal heights 
measured directly from the 3D model snapshot. The tem- 
perature variations {5T) are defined relative to the aver- 
age temperature (T) over the snapshot at each particular 
height, i.e. ST = {T~f)/f. It is clear from the figure that 
the dependences measured directly from the 3D model are 
very similar to the ones obtained after spectral synthesis 
and the A-meter technique (see dashed-dotted curved with 
open circles on the same figure). The only remarkable dif- 
ference is that the A-meter velocities and contrast values 
are systematically lower. The latter is easy to understand 
having in mind that the A-meter method gives the aver- 
age information over a certain height range, thus leading to 
such a reduction. 

6. Discussion 

The classical picture of the convection breaks in the solar 
photosphere, starting at least from 50 km above the contin- 
uum formation level. The convective elements moving up- 



wards and downwards not only change their sign of contrast 
but also their direction of motion. Often both events take 
place at the same time. This process occurs in a wide range 
of heights up to 650 km. The height where the contrast 
sign reversal takes places depends strongly on the contrast 
of gr anules arid inte rgranules at the continuum formation 
level. IStodilkal (|2006D came to the conclusion that the tem- 
perature inversion happens in a wide range of heights, from 
150 to 500 km. In addition, there is a similar tendency, i.e. 
the larger the granule, the h igher the temperature i nver- 
sion takes place. According to lPuschmann et all (|2005[ ). the 
temperature inversion occurs for structures with sizes over 
l."5 at heights above 200 k m. Our results co nfirm these con- 
clusions. On the other hand. lStodilkal (|2006D claims that the 
temperature reversal events are more often above intergran- 
ules, while according to our results, this happen more often 
above granules (see right panel of Fig. [9]) . The large scatter 
in the average height of the contrast (temperature) reversal 
according to observations of different authors (from 60 to 
350 km) not only stems from the observational details (like 
seeing) and their analysis, but also from the methods used 
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out the contrast sign reve rsal i n the solar atmosphere wer e 
lEvans fc Catalanol (|1972D and iHolweger fc Kneerj ()1989[ ). 
However, not only does the contrast sign reversal take place 
in the solar atmosphere, but also the reversal of sign of the 
convective velocities with height. On average, it happens at 
about 200-300 km high. As far as we know, this property 
of the convective motions has not ever been studied before. 

The comparison of observed dependences with those ob- 
tained from the NLTE synthesis of Bail 4554 A profiles 
in the 3D hydrodynamical model of lAsplund et al.l (|2000r ) 
have shown that this model reproduces the observed depen- 
dences rather well except for the dependence of the height 
of the contrast sign reversal on the velocity at the bottom 
layer. We found that all 16 classes of the convective motions 
are present in the 3D model, and their statistical behavior is 
similar to the observed. The contrast sign reversal and the 
velocity sign reversal take place over the whole atmosphere 
from 50 to 650 km. However, the column structure of gran- 
ulation is preserved up to 650 km. Almost 40% of all the 
convective structures detected at the continuum formation 
height reach 650 km. In particular, 80% of the structures 
larger than 1500 km in the continuum reach 650 km. Most 
losses are experienced by the smallest convective elements 
with sizes below 500 km, while the fraction of the elements 
with larger sizes increases with height. 



to calculate the formation height s of the observed spectra l 
lines, as correctly pointed out in lPuschmann et all (|2005l ). 
We think that the main reason for this disagreement is the 
difference in the formation heights of spectral lines used in 
different observations. To prove this statement, consider the 
calculation presented in Fig. [T31 The panel a of this figure 
gives the average height of the contrast sign reversal as a 
function of the Ban line depth for the our observations. It 
follows that the larger the line depth, the higher the level of 
the contrast sign reversal obtained. The panel b of Fig. [13] 
gives the fraction of the contrast sign reversal events as 
a function of the line depth. As one can see, the stronger 
lines probe higher atmospheric regions and may, therefore, 
detect more sign reversals. Since different authors use spec- 
tral lines of different strengths (corresponding to different 
formation heights of the central intensity), it is natural 
that they obtain different heights of the sign reversal of 
solar granulation contrast. The panels c and d of Fig. [T3] 
give similar results but for the heights where the velocity 
sign reversal takes place. The first authors that pointed 



7. Conclusions 

We performed statistical analysis of the solar granulation 
velocities and intensities extracted by means of a A-meter 
technique from the high-resolution observations of the quiet 
solar disk center in the Ban 4554 A line. As this line forms 
high enough in the atmosphere, it has allowed us to study 
the properties of granulation at heights not accessed in the 
earlier studies up to 650 km, i.e. the lower chromosphere. 
The results of observations were compared with the syn- 
thetic data fro m the 3D hydro dynamical model of the solar 
convection bv lAsplund et al] f2000). We found similar de- 
pendencies in the observations and simulations smeared to 
the observational resolution. It gives us confidence in our 
results and also proves the good quality of the 3D model. 
Our main findings can be summarized as follows 

— About 40% of the column structure of granulation is 
preserved up to heights around 650 km. 
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— The average size of granular elements increases with 
height. 

— The granulation can be described by a 16-column model, 

depending on the sign of contrast and velocity direction 
at the bottom and top atmospheric levels. All types of 
motion are found to indeed be present on the Sun. 

— The contrast sign reversal of granulation takes place at 
heights around 200-300 km. At the same heights, on 
average, the velocity sign reversal also occurs. 

— The heights where the sign reversal happens strongly 
depends on the individual contrast and velocity of gran- 
ules at the continuum. The larger these parameters, the 
higher the sign reversal takes place. 

— The precise values of the heights of velocity and contrast 
sign reversal obtained from observations are sensitive to 
the spatial resolution within their standard deviation 
error bars. 
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